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ABSTRACT 

(Distribution  Limitation  Statement  B) 


The  heating  of  a  solid  body  by  a  laser  beam  has  been  analyzed  using  a  generalized 
two-dimensional  solution  of  the  heat  conduction  equation  for  a  finite  right 
circular  cylinder.  This  series-eigenvalue  solution  Includes  time-dependent  flux 
boundary  conditions  and  reradiation  from  all  surfaces.  Temperature  solutlins  are 
derived  for  three  spaclal  laser  beam  profiles  applied  to  continuous  wave*  pulsed* 
and  repetitively  pulsed  laser  heating.  Computer  programs  developed  to  evaluate 
these  temperature  solutions  are  presented. 
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SECTION  I 
INTRODUCTION 

The  Investigation  of  the  effects  of  continuous,  pulsed,  and  repetitively 
pulsed  lasers  on  materials  has  brought  about  a  renewed  effort  to  solve  a  partic¬ 
ular  class  of  heat  flow  problems.  These  problems  vary  greatly  In  complexity 
depending  on  the  thermal  properties  of  the  material  of  Interest;  the  temporal 
history  of  the  laser  flux;  the  spaclal  Intensity  profile  of  the  laser  beam;  and 
the  shape,  size,  and  optical  properties  of  the  body  being  Irradiated.  Few  laser 
heating  problems  can  be  considered  steady  state,  and  most  are  two-  or  three- 
dimensional  In  nature. 

Ready  (Ref.  1)  has  compiled  an  excellent  summary  of  heat  transfer  solutions 
applicable  primarily  to  pulsed  laser  heating  problems.  There  are  solutions  for 
two-dimensional  radially  symmetric  heating  by  circular  and  Gaussian  beam  pro¬ 
files  as  well  as  one-dimensional  cases.  Most  of  these  have  been  solved  for  the 
semi-infinite  solid  case  because  of  the  emphasis  on  pulsed  laser  heating, 
although  some  solutions  are  presented  for  the  finite  thickness  plate.  Including 
surface  losses  from  reradiation  or  convection.  To  obtain  these  concise  analyt¬ 
ical  solutions,  It  was  assumed  that  the  heat  flow  obeys  the  classical  heat 
transfer  equations  with  no  phase  changes  occurring.  Material  thermal  properties 
are  Independent  of  temperature.  Despite  further  assumptions  about  the  temporal 
shape  of  the  laser  pulse  and  small  absorption  depths  at  surface  (excellent 
assumption  for  metals),  some  of  the  solutions  presented  still  require  the 
numerical  evaluation  of  Integrals. 

Carslaw  and  Jaeger  (Ref.  2)  give  a  one-dimensional  transient  heat  transfer 
solution  for  a  finite  thickness  plate  with  a  prescribed  flux  entering  one  surface 
and  no  heat  flow  across  the  other  surface.  Again,  constant  thermal  properties 
were  assumed,  and  no  phase  changes  are  permitted.  If  absorption  lengths  for  the 
laser  energy  are  assumed  to  be  short,  this  solution  can  be  used  to  model  contin¬ 
uous  wave  laser  heating  with  the  condition  that  the  one-dimensional  criteria  are 
met.  This  requires  that  for  a  material  of  given  thermal  properties,  the  average 
laser  power  be  sufficiently  high  and  the  sample  be  sufficiently  thin  (Ref.  3). 
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An  additional  one-dimensional  heat  transfer  solution  for  a  finite  thickness 
plate  applicable  to  pulsed  laser  heating  can  be  derived  by  (1)  using  the  above 
prescribed  flux  solution  for  the  duration  of  the  pulse,  and  (2)  applying  the 
resulting  temperature  profile  at  the  end  of  the  pulse  as  an  initial  condition  to 
a  finite  thickness  plate  solution  (Ref.  2)  for  any  Initial  temperature  profile 
and  both  surfaces  Insulated.  This  second  solution  then  describes  the  diffusion 
away  from  the  front  surface  of  the  heat  deposited  by  the  pulse.  Again  the  same 
criteria  for  one-dimensionality  must  be  obeyed. 

Since  In  most  laser  experiments  the  one-dimensional  requirements  are  not  met. 
It  Is  desirable  to  have  a  full  transient  two-dimensional  heat  transfer  solution 
to  theoretically  describe  laser  heating.  N.  Y.  Olcer  (Refs.  4,5)  has  derived 
very  general  solutions  of  the  appropriate  heat  conduction  equation  for  the  heat¬ 
ing  of  a  finite  solid  body  of  arbitrary  shape,  any  Initial  temperature  profile, 
and  time  dependent  surface  conditions.  In  this  work  he  has  specialized  the 
solution  to  that  for  a  finite  right  circular  cylinder  that  retains  the  very 
general  time  dependent  boundary  conditions  which  permit  his  results  to  be 
applied  with  appropriate  simplification  to  the  laser  heating  problem. 

It  Is  the  Intent  of  this  report  to  show  how  the  work  of  Olcer  can  be  simpli¬ 
fied  and  yet  be  made  to  model  laser  heating  In  as  general  a  way  as  possible. 
Specific  examples  are  then  derived  for  heating  with  three  different  spaclal 
laser  beam  Intensity  profiles  and  for  continuous  wave,  pulsed,  and  repetitively 
pulsed  lasers.  The  computer  programs  used  to  evaluate  the  solutions  are 
described,  and  two  examples  are  given  In  the  appendix. 

Finally,  the  applications  section  outlines  ways  to  use  the  existing  programmed 
solutions.  It  should  be  remembered  that  many  other  solutions  for  laser  heating 
problems  can  be  derived  from  Olcer's  basic  solution  In  addition  to  the  ones 
presented  here. 
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SECTION  II 

APPLICATION  OF  OLCER'S  GENERAL  SOLUTION  TO  THE  LASER  HEATING  PROBLEM 

* 

1.  SIMPLIFICATION  OF  OLCER'S  SOLUTION 
The  heat  conduction  equation 


V’T'Or'.t')  +  Q'dr'.t') 


1 

K 


3T'(r',t') 

at1 


r'  In  R' 
t'  >  0 


0) 


for  the  heating  of  a  finite  Isotropic,  homogeneous  solid  body,  R' ,  of  arbitrary 

•• 

shape  with  constant  thermal  properties  has  been  solved  by  N.  Y.  Olcer  (Refs.  4, 
5).  The  following  general  boundary  conditions  apply  on  the  surfaces 
SJ(1*l,2,...,q)  of  the  body  R', 


t'  >  0 

r'  on  Sj(1*l,2,...,q) 


(2) 


with  an  Initial  temperature  distribution 


T* (r * ,t ' )  *  F'(r')  r'  In  R'  and  on  Sj ;  t'  s  0  (3) 

The  terms  of  equations  (1),  (2),  and  (3)  are  defined  as  follows: 

T'(r',t')  *  unsteady  temperature  distribution  In  the  region  R' 

V'2  *  the  Laplaclan  in  r'-space 

Q'(r',t')  ■  heat  production  per  unit  time  per  unit  volume  In  R' 

K  *  thermal  conductivity 
k  •  thermal  dlffuslvlty 
r*  *  position  of  a  point  in  the  region  R' 
t'  *  time 

q  *  number  of  surfaces  Sj  bounding  R' 

fj(r' ,t')  »  heat  entering  surfaces  Sj  per  unit  time  per 
unit  area 
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hj(r')  *  a  coefficient  of  heat  transfer  (radiation  and 
convection)  from  surfaces 

n^  *  direction  of  the  outward  normal  to  surface 

F'(r')  *  prescribed  Initial  temperature  distribution 

For  the  convenience  of  removing  constants  from  equations  (1)  and  (2)  and 
simplifying  Its  solution,  the  following  nondlmenslonal  parameters  are  defined: 

ft(r,t)  *  a  fj(r',t')/T0K 

h^r)  ■  a  hj(r')/K 

Q(r,t)  *  a2Q' (r'.fj/y 

r  a  r'/a 

t  =  fic/a2 

T(r,t)  -  T'(r\t,)/T0 

F(r)  -  F'(r')/T0  (4) 

where  a  Is  a  characteristic  dimension  of  the  body  R'  and  TQ  Is  some  reference 
temperature.  Equations  (1),  (2),  and  (3)  can  now  be  written  In  dimensionless 
form. 


V2  T(r,t)  +  Q(r,t)  - 

3~+  (r)l  T(r,t)  =  f^r.t) 

T(r,t)  -  F(r) 


r  In  R 
t  >  0 


(5) 


r  on  S^(1-l,2,...,q) 
t  >  0 


r  In  R  and  on 
t  -  0 


(7) 


t* 
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where  R,  S^,  n^,  and  V2  In  nondlmenslonal  r-space  correspond  to  th®  same 
quantities  or  operations  primed  In  r'-space  with  dimensions.  All  further 
discussion  of  solutions  of  the  heat  transfer  equation  In  this  report  will  be 
In  terms  of  dimensionless,  unprimed  quantities.  In  the  numerical  evaluation 
of  a  solution  equation  (4)  must  be  used  first  to  translate  Input  parameters 
to  dimensionless  form,  and  second  to  return  the  calculated  temperature  field 
to  actual  temperatures. 

The  solution  to  equations  (5),  (6),  and  (7)  Is  from  reference  5. 


T(r,t)  *  £  Am  *m(r)exp(-y£t) 
m 


///  «’m(?)F(7)dR 


+  f  exp(y2x)  ffm 

ff  ^Jr)Q(rt  r)dR  + 

&  //  *n<?>fi<?-T>dsiV 

0  \  F 

i 

n  S1  /  1 

where  the  eigenfunctions  ^m(r)  and  the  non-negative  eigenvalues  y^  are  o 
by  solving  the  eigenvalue  problem 

(V2  *  vm2 

)  4»m(r)  *  0 

r  In  R 

feth*rr! 

)  *m(r>  =  0 

r  on  S^(1-l,2,...,q) 

(9) 


and 


III 


r  ■  JiJ  *m(7,dR 

m  n 


(10) 


This  solution  has  been  specialized  to  the  case  of  a  right  circular  cylinder 
of  radius,  a,  and  thickness,  2b,  In  cylindrical  coordinates  (Ref.  5).  For  the 
application  of  this  solution  to  the  heating  of  a  disc  by  a  laser.  It  will  be 
assumed  that  the  spaclal  Intensity  profile  of  the  laser  beam  Is  symmetrical 
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about  Its  axis,  and  therefore,  the  temperature  profile  In  the  heated  body  also 
has  cylindrical  symmetry.  Dropping  the  angular  dependence  the  nondlmenslonal 
coordinates  required  are 


and 


r  -  r'/a;  z  -  z'/a 


n  =  b/a 


The  heat  flow  equation  to  be  solved  and  the  necessary  boundary  and  Initial 
conditions  corresponding  to  equations  (5),  (6),  and  (7)  are  as  follows. 


— JT(r,z,t)  +  Q(r,z,t)  *  0 
at/ 

0  <  r  <  1 


|z|  <  n  t  >  0  (11) 


(“  Iz  +  hi)  T(r»z»t)  * 

Oi  r<l  z«-n  t>0 

(fz  +  hz)  T(r»z»t)  *  f2(r.t) 

Oir  < 1  z*n  t>0 

(iz  +  hs)  T^r»z*t)  *  Mz»t) 

r  ■  1  |z|  <  n  t  >  0  (12) 


T(r,z,t)  *  F(r,z) 

0<r<l  |z|  <  n  t-0  (13) 

The  solution  of  equations  (11),  (12),  and  (13)  for  the  axlsymmetrlc  right 
circular  cylinder  which  corresponds  to  equation  (8)  Is 


6 


T(r,z,t)  =  ti  cl.Lk*n<r)xk(z)  e*P  [-(x^k)‘] 


k  n 
L  n 


<j>n(r)xk(z)F(r,z)rdrdz  + 


/  "p  [(x^k)T] 


o  rn 
■  i  n 


i 

*n(r)xk(z)Q(r,z,T)rdrdz  +  j  <*n(r)  (f^v)  +  xk(o)f2(r,t))  rdr 


f  xk(z)f3(z,t)dz  dt 


and  the  eigenvalue  problem  of  equation  (9)  becomes 


♦  II.* 

r  3r  3z2  j 

)  V>z)  =  0 

Oir  <  1 

N 

A 

J 

(-l?+h.) 

|  *kn(r.z)  ■  0 

0ir<  1 

z  ■  -n 

(l?  *  % 

|  *kn(r.z)  -  0 

Oir  <  1 

z  *  n 

(K 

)  'Pkn(r-Z)  ’  0 

r  *  1 

|z|  <  n 

(15) 

The  eigenfunctions  ^n(r,z)  and  the  eigenvalues  ukn  are  double  Index  sets  where 
k,n  =  0,1,2,...“.  Where  hj.hj.h,  are  constants,  the  solution  to  equation  (15) 
Is  a  product  of  two  eigenfunctions 

*k,n(r>z)  '  *n(r)*k(z)  061 

where  the  eigenfunction  <j>n(r)  Is 

*„(D  -  J,  (V)/J.(X")  071 
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(18) 


and  X  is  the  nth  non-negative  root  of 

Mo(Xn)  '  xn  Ji(xn) 

The  eigenfunction  XjJz)  is 


and  vk  Is  the  kth  non-negative  root  of 

(hj+h2)vk  cos  2nvk  =  (v£-h,h2)  s1n  2rivk  (20) 

and  the  two  eigenvalue  sets  are  related  by 


Mkn  <"> 

From  equation  (10) 


B  2r)  for  vk  *  0  (23) 


The  temperature  expression  of  equation  (14)  and  Its  attendant  definitions 
(equations  (17),  (18),  (19),  (20),  (22),  and  (23))  are  still  much  more  general 
than  necessary  to  solve  the  laser  heating  problem.  A  few  observations  will 
reduce  equation  (14)  to  a  much  more  manageable  size.  First,  In  the  typical 
laser  experiment  there  Is  Initially  no  temperature  profile  In  the  body  being 
Irradiated,  so  that  Integral  containing  the  function  F(r,z)  can  be  dropped. 


8 


Second,  by  assuming  that  the  material  is  a  good  absorber  of  laser  energy,  we 
are  Irradiating  only  the  front  surface  of  the  cylinder  z  =  n  with  flux  f2(r,t), 
and  therefore,  fj(r,t)  and  f3(z,t)  are  both  zero.  The  Integral  containing 
Q(r,z,r)  also  vanishes  because  no  heat  Is  being  generated  in  the  bulk.*  Third, 

i 

by  assuming  that  the  coefficients  of  combined  reradiation  and  convection  are 
the  same  for  all  cylinder  surfaces  h1=h2*hJ*h,  equations  (20)  and  (23)  are 
simplified  slightly.  We  now  rewrite  equation  (14)  with  Its  accompanying 
definitions  for  the  laser  heating  problem. 

00  00 

T(r,2‘t)  '  ?  ?  CnLkVr>X|c(z)  exp[- (x>k)  ‘] 


}  r  r 

1 

J  J  exe  [(^k)1 
0  ^0  L 

r  d>n(r)  xk(n)f2(r,T)rdrdt 

where 

*„<''>  ■  J.(vVJ.(x")  (25) 

An  is  the  nth  non-negative  root  of 

*  \  (26) 


xk(z)  =  cos  jvk(n+z)J  +  (h/vk)s1n  jvk(n+z)l 
vk  Is  the  kth  non-negative  root  of 

2hv^cos  2nvj(  =  ^-h2^sin  2nvjc 


(27) 


(28) 


♦For  dielectric  materials  with  very  long  absorption  lengths  at  the  laser  wave- 
lenth  of  Interest,  the  absorption  of  energv  In  the  bulk  could  be  handled  by 
retaining  the  heat  generation  term  Q(r,z,T)  and  setting  the  three  surface 
fluxes  fj(r,t),  f2(r,x),  and  fs(z,t)  to  zero. 
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Cn  *  MXn+h!)  («> 

Lk  =  vk/[,1(vk+h2)+  h]  vk  *  0 

■  l/2n  vk  “  0  (30) 

The  temperature  expression  summarized  by  equations  (24)  through  (30)  Is 
still  a  very  general  solution  for  the  heating  of  a  cylinder  by  a  laser,  because 
the  only  requirement  made  Is  that  the  flux  entering  the  front  surface  be  radially 
symmetric.  The  flux,  represented  by  f2(r,i),  Is  still  a  general  function  of 
radius  and  time.  It  Is  desirable,  of  course,  to  choose  beam  Intensity  profiles 
In  time  and  space  which  when  substituted  Into  equation  (24)  permit  a  direct 
analytical  evaluation  of  the  double  Integral  In  the  braces  {  }  because  this 
Integration  Is  Inside  of  the  double  sum  and  must  be  performed  n»k  times. 

2.  LASER  BEAM  SPACIAL  INTENSITY  PROFILES 

The  spaclal  Integration  of  equation  (24)  has  been  evaluated  for  three  laser 
beam  profiles:  constant  Intensity,  fj,  of  finite  radius  (cylindrical  profile), 
Gaussian,  and  parabolic.  For  the  cylindrical  profile, 

f2(r)  =  f;  0  <  r'  <  o' 

f2(r)  =0  o'  <  r' 

and  the  total  beam  power,  P',  Is 

P'=TO,2f;  (31) 

Normalizing  the  appropriate  parameters  the  Integral  to  be  performed  Is 
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For  the  Gaussian  profile 

fj(r)  =  fj  exp(-r'2/2a'2)  0  <.  r*  ia1 

and 


P'  =  27Ta,2f: 


(33) 


With  normalized  parameters  the  required  integral  is 


/ 


exp(-r2/2o2)rdr 


This  Integration  cannot  be  completed  In  closed  form;  however,  an  exact  expres¬ 
sion  is  available  for  the  definite  integral  from  zero  to  infinity. 


oo 

/ 


J  ,(6r)rv+1  exp(-y2r2)dr  * 

V  (2y2) 


6' 


exp(-B2/4y2) 


(34) 


As  long  as  the  effective  radius  of  the  beam  is  somewhat  smaller  than  the  sample 
radius,  the  exponential  part  of  the  function  to  be  Integrated  has  decayed  to 
negligible  values  for  r  *  1,  and  thus,  integrating  to  1  is  equivalent  to 
integrating  to  Applying  equation  (34)  for  v  =  0, 


i 


exp(-r2/2a2)rdr  ^ 


(35) 


For  the  parabolic  profile 

fj(r)  =  fj(l  -  r,2/0*2)  0  <  r'  <  o' 

fj(r)  =0  o'  <  r' 


and 


P'  =  7ra,2f;/2 


(36) 


11 


ArwL-iK-/z-z;tt> 


The  required  Integral  Is  evaluated  with  an  integration  by  parts 


c  i 

J,  Jo 


W  ' 


(1  -  r2/a2)rdr 


2f.  J, 
0  2 


M 


^  J.(Xn) 


(37) 


3.  LASER  BEAM  TEMPORAL  INTENSITY  PROFILES 

The  temporal  integration  of  equation  (24)  has  been  performed  for  three 
different  kinds  of  lasers,  continuous  wave  (CW),  pulsed,  and  repetitively 
pulsed  (REP).  The  REP  solution  works  for  a  very  general  pulse  shape  as  each 
pulse  Is  represented  by  a  rising  and  a  decaying  exponential.  Two  types  of 
single  pulse  solutions  are  easily  derived.  The  first,  a  simple  square  pulse 
Is  a  specialization  of  the  CW  solution,  while  the  second  Is  obtained  directly 
from  the  REP  solution  by  setting  the  pulse  frequency  to  a  number  sufficiently 
small  to  allow  only  one  pulse  of  energy  during  the  times  of  Interest. 

For  the  CW  laser  beam 


fi(t')  -  f,  0  <  t' 

and  the  time  integration  in  equation  (24)  is  just 

j[  exp  [fatuk)T]  f.dT  "  77^7  jexp  [ftS)1]  - 1 1  (38) 

The  single  square  pulse  beam 

fj(t')  ■  f0  0  <  t'  <  tp 

f;(f)-0  t'<f 

requires  only  a  change  of  the  upper  limit  of  equation  (38)  to,  tp,  the  normal¬ 
ized  pulse  length.  (38a) 

A  temperature  solution  Is  possible  for  the  REP  laser  If  It  Is  assumed  that 
the  pulses  are  all  of  Identical  energy,  shape,  and  spacing  In  time.  The  follow¬ 
ing  general  shape  was  chosen  for  a  pulse. 
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0  <  t'  <  t' 

-  upp 


rise:  f;(t')  «  fj  [l -exp  (-t7t'pr)] 

decay:  fj(t')  =  fj  exp[-(t'-fpp)/t'pd]  ‘pp  ' 

Three  time  constants  are  required  to  describe  each  pulse:  the  rise  time,  tpr, 
the  time,  tpp,  for  the  Intensity  to  reach  a  maximum,  fj,  and  the  decay  time, 
tpd*  Tfie  ener9y  Per  area  contained  In  this  pulse  Is  given  by 

nE?ir9y — =  f '  ft*  -t'  +t'j+t'  exp/-t'  /t'  VI  (39) 

Unit  area  pp  pr  pd  ‘'pr  pp'  pr/J  v  ' 

Thus,  for  the  mth  pulse  In  a  series  of  pulses  separated  by  a  period,  81, 
rise:  f ^ ( t ' )  *  fj[l-exp(- (t'-(m-l)6')/t^r)J 

(m-l)e'  <  t'  £  (m-l)e'  +  t^ 

decay:  fj(t')  -  f;  exp[- (t*-(m-l Je’-t^/t^] 

(m-l)e'  +  t^p  <  t'  i  me' 

The  time  Integral  of  equation  (24)  where  t'  Is  a  time  during  the  Mth  pulse 

t 

/  exP[(Xn+vk)T]f2(T>dT 

0 

therefore,  consists  of  first,  a  complete  Integration  over  pulses  1  through  M-l, 
and  second.  Integration  over  the  Mth  pulse  up  to  time  t'.  Substituting  a 
normalized  form  of  f^  (t ' )  for  pulses  1  through  M-l  and  for  pulse  M  Into  the 
time  Integral  and  letting  the  sum  of  the  squares  of  the  elgenvlaues,  A*+v£, 
be  represented  by  Enk  one  obtains 
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F^rr  hP  [(Enk'Vtpr)tpp]  ‘  \ 


either 


or 


+  E^TT  “*ta/**)[«*  IMS*] 

-  e*p[(Enk-1/tpd)tpp]]|  0  <  t  i  (M-l)6 

+  rJik  '  exp  [E"k  ((M-l  )e-t)]  J 

'  [exp  )6‘t)/tpr]  -  '*p[e„I<(<h-’H]] 

(M-i)e  <  t  <  (M-i )e  +  t 
(a  time  during  the  rise  of  the  Mth  pulse) 

+  exP[Enk (<M-1  >«-»)]  |e^  [«p(Enktpp)-l] 

■  h  ^k-’^pr)^]-’]  ( 

+  I^k^dIT  [exP  [((M-,)e+tpp-t)/tpd]-  «p[E„k((«-i)enpp-t)]] 


(M-i)e  +  tpp  <  t  <  Me 

(a  time  during  the  decay  of  the  Mth  pulse) 

(40) 
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The  first  part  of  equation  (40)  containing  ^  Is  easier  to  evaluate  than 
It  appears  at  first  because  the  Index,  m,  and  the  variable,  t,  are  contained 
only  In  the  first  exponential,  term,  and  the  argument  Is  always  negative.  Since 
REP  lasers  to  date  operate  In  the  order  of  hundreds  of  pulses  per  second  or 
less,  M,  the  total  number  of  pulses  Is  never  an  unreasonably  large  number. 
Furthermore  the  expression  contained  In  the  braces  Is  just  a  double  Indexed 
set  of  constants  to  be  evaluated  once  after  the  eigenvalues  are  found.  The 
choice  of  the  appropriate  second  part  of  equation  (40)  Is  made  directly  from 
the  variable  t. 

As  stated  earlier,  this  solution  can  be  used  In  the  single  pulse  case  by 
entering  a  period,  0,  which  Is  always  greater  than  t.  Only  the  appropriate 
second  part  of  equation  (40)  Is  to  be  evaluated  for  M*l. 
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SECTION  III 
THE  PROGRAM 

To  obtain  a  solution  to  equation  (24)  one  must  consider  sufficiently  many 
terms  of  the  series  to  allow  for  convergence.  In  general,  the  number  of  terms 
required  for  convergence  Is  large  enough  to  make  any  solution  other  than  by  a 
high-speed  computer  Impractical.  The  program  which  Is  described  below  and 
appears  In  the  appendix  was  written  to  provide  an  evaluation  of  equation  (24). 
This  program  was  written  In  FORTRAN  IV  source  language  and  Is  run  on  a  CDC 
6600.  The  following  Is  a  list  of  symbols  used  In  the  program. 


CAP 

C0ND 

A 

D 

SIG 

FL 

H 

ELAM 

ENU 

P 

TPR 

TPP 

TPD 

EGN 


K 

K 

a 

2b 

o' 

n 

h* 

X 


v 


P' 


*PP 

lpd 


The  main  program,  called  TMPC,  Initializes  Input  data  through  the  use  of 
data  statements.  The  desired  times  of  Interest  are  generated  prior  to  any  call 
to  subroutines.  In  addition,  one  takes  Into  consideration  the  relationship 


u 


iv 


1 
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between  peak  intensity  and  power  as  described  in  equations  (31),  (33),  and  (36) 
according  to  the  particular  spaclal  and  temporal  beam  profile  of  interest.* 

After  these  preliminary  values  have  been  determined,  then  a  call  Is  made  to 
INITL.  This  subroutine  first  computes  the  heat  transfer  coefficient  H,  and  then 
normalizes  various  constants.  Finally,  two  calls  to  the  subroutine  ZER0  are 
made:  one  for  the  function  TRANT,  and  one  for  the  function  TRANB. 

Subroutine  ZER0  determines  a  preassigned  number  of  roots  to  a  desired 
accuracy  for  a  given  function.  The  function  TRANT  evaluates  equation  (28)  for 
the  particular  value  of  v  which  Is  passed  to  It.  Slmllarily,  TRANB  evaluates 
equation  (26)  for  different  values  of  X.  The  number  of  v^'s  and  Xn's  calculated 
Is  determined  by  NEGT  and  NEGB,  respectively.  It  has  been  found  that  2000  vk's 
and  300  An's  will  normally  produce  a  temperature  accuracy  within  3  x  10“** 
degrees  centigrade  with  the  time  required  to  calculate  them,  usually  In  the 
range  of  a  fraction  of  a  second.  If  computer  time  becomes  excessive  these 
numbers.  In  many  cases,  can  be  substantially  decreased.  The  number  of  terms 
actually  needed  Is  determined  primarily  by  the  magnitude  of  the  temperature 
gradient  at  the  particular  time  of  Interest.  For  example,  the  number  of  Xn's 
required  for  the  cylindrical  beam  at  short  times  Is  significantly  greater  than 
the  number  required  In  the  Gaussian  case  at  the  same  times.  After  the  vk's 
and  the  Xn's  have  been  calculated  control  Is  returned  to  the  main  program. 

A  call  Is  then  made  to  the  function  routine  TMP  which  will  compute  the 
temperature,  as  given  In  equation  (24),  for  a  particular  radius,  depth,  and  time. 
In  TMP  a  desired  accuracy  Is  set  for  the  temperature  calculation,  usually 
3  x  10"4,  and  the  radius,  depth,  and  time  are  normalized.  The  function  then 
begins  to  evaluate  the  first  term  of  the  double  sum  of  equation  (24)  by  setting 
the  Indices  k  and  n  to  1.  ARGA  and  ARGB  are  then  calculated  to  be  later  used 
as  the  arguments  In  the  computation  of  $  as  defined  In  equation  (25).  Next 
a  call  Is  made  to  function  routine  SINT  which  does  the  spaclal  Integration  of 
equation  (24)  for  either  a  cylindrical  profile  (equation  (32)),  Gaussian  profile 
(equation  (35)),  or  parabolic  profile  (equation  (37)).  The  program  next  com¬ 
putes  the  k,nth  term  of  the  double  sum  In  the  following  order. 


♦There  are  two  sample  programs  listed  In  the  appendix:  one  for  a  CW  laser  with 
a  cylindrical  beam  profile,  and  another  for  a  REP  laser  with  a  Gaussian  beam 
profile. 
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First,  the  functional  constant  Cn,  evaluated  In  function  routine  C,  is 
determined  as  defined  In  equation  (29).  The  next  term  Is  the  functional  con¬ 
stant  Lk,  evaluated  In  function  routine  SL,  as  defined  In  equation  (30).  The 
previously  determined  arguments,  ARGA  and  ARGB,  are  used  In  the  computation  of 
the  eigenfunction,  4>n.  The  following  two  terms  of  the  product  are  functions 
found  In  equation  (26)  and  evaluated  In  the  function  routine  CHI. 

The  product  of  the  temporal  Integration  and  the  time  dependent  exponential 
term  of  equation  (24)  computed  In  function  routine  TFACT,  comprises  the  next 
term  of  the  product.  The  temporal  Integration  will  be  given  by  either  equation 
(38),  (38a),  or  (40)  according  to  whether  the  laser  beam  Is  CW,  single  pulse, 
or  multiple  pulse.  Here  "a"  denotes  a  change  In  the  upper  limit  of  equation 
(38)  to  tp.  The  final  factor  Is  the  result  of  the  spaclal  Integration  pre¬ 
viously  computed.  The  program  will  continue  to  run  through  this  double  sum 
until  either  the  specified  accuracy  Is  achieved  or  the  number  of  eigenvalues 
requested  Is  exceeded.  Once  one  of  those  conditions  has  been  met.  In  the  latter 
case  an  error  message  will  be  printed,  then  control  Is  returned  to  the  main 
program.  Repeated  calls  to  TMP  will  produce  the  desired  temperatures  and  all 
that  remains  Is  the  output  of  this  Information. 

A  typical  output  from  this  program  might  Include  thermal  properties  of  the 
substance  under  consideration,  dimensions  of  the  sample,  beam  parameters,  and 
the  respective  times  and  temperatures  at  specified  locations. 

The  program  usually  will  compute  10  temperatures  on  a  CDC  6600  In  about 
1  or  2  seconds. 
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SECTION  IV 
APPLICATIONS 

1 .  DIRECT  CALCULATION  OF  TEMPERATURES 

Direct  calculation  of  temperature  time  histories  at  any  number  of  locations 
or  temperature  profiles  at  any  times  are  possible  with  a  minimum  expenditure  of 
computer  time.  Because  this  Is  a  series  analytical  solution  it  Is  necessary  to 
calculate  only  the  temperatures  at  the  locations  and  times  desired,  rather  than 
generate  the  whole  temperature  field  at  every  time  step  up  to  the  time  of 
Interest  as  done  In  finite  element  or  finite  difference  computer  solutions. 

This  analytical  solution  Is  particularly  useful  for  calculating  heat  flow  In  a 
variety  of  laser  effects  experiments.  Temperature  rises  can  be  estimated  for 
two-dimensional  CW,  pulsed,  and  REP  lasers.  The  effect  of  temporal  pulse  shape 
on  the  temperature  history  can  be  Investigated  In  detail,  and  In  the  case  of 
REP  lasers,  the  amount  of  radial  and  axial  "recovery  cooling"  of  the  Irradiated 
area  occurring  between  pulses  can  also  be  calculated. 

In  addition  to  the  Investigation  of  various  time  parameters  characterizing 
the  laser  flux,  the  Importance  of  the  laser  beam's  spaclal  profile  can  also  be 
studied.  For  a  laser  beam  of  a  given  total  power,  the  beam  size  parameter,  a, 
and  peak  Intensity,  f0,  can  be  varied  to  determine  the  effects  of  focusing  or 
defocuslng  the  laser  beam  on  the  temperature  profiles  In  a  sample.  Also  tem¬ 
perature  fields  created  by  the  different  beam  spaclal  profiles  can  be  compared, 
at  least  for  the  cylindrical,  Gaussian,  and  parabolic  profiles  derived  here. 
Certainly  other  mathematical  functions  can  be  chosen  to  represent  the  beam 
profile  of  Interest.  With  several  possible  functions  to  try,  one  can  compare 
the  resulting  temperature  field  and  determine  quantitatively  how  Important  Is 
the  precise  representation  of  che  laser  beam  profile  to  the  predicted  form  of 
the  temperature  field. 

2.  DETERMINATION  OF  SURFACE  ABSORPTION  COEFFICIENTS 

In  the  study  of  the  effects  of  laser  beams  on  materials  It  Is  nearly  always 
desirable  to  know  what  fraction  of  the  Incoming  laser  beam  Is  being  absorbed  by 
the  sample.  This  surface  absorption  coefficient  can  be  obtained  from  experi¬ 
mental  temperature  time  histories  with  the  aid  of  this  heat  transfer  solution. 
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For  example,  temperature  data  from  thermocouples  attached  to  a  laser  Irradiated 
sample  can  be  fit  by  a  method  of  "least  squares"  with  the  mathematical  form  of 
the  temperature  rise  predicted  by  the  two-dimensional  program.  Knowing  the 
actual  beam  power,  the  program  can  output  what  fraction  of  that  power  the 
sample  must  have  absorbed. 

This  data  fitting  technique  permits  one  to  obtain  surface  absorption  coeffi¬ 
cients  for  the  wide  range  of  laser  effects  experiments  where  the  assumption  of 
one-dimensional  heat  flow  cannot  be  made.  In  fact  the  amount  of  radial  heat 
losses  In  most  laser  experiments  Is  considerable,  and  If  not  accounted  for 
properly  can  result  In  enormous  errors  In  experimental  absorption  coefficients. 

3.  OTHER  APPLICATIONS 

Because  of  the  remarkably  general  nature  of  this  two-dimensional  heat 
transfer  solution  as  expressed  In  equation  (14),  one  can  envision  a  variety  of 
other  applications.  Rather  than  assuming  surface  heating,  the  laser  energy 
could  be  absorbed  over  finite  depths  In  the  sample  through  the  use  of  the  heat 
generation  Integral.  More  sophisticated  solutions  might  then  be  developed 
which  would  model  the  heating  of  materials  semitransparent  to  laser  radiation 
such  as  lenses,  windows,  etc. 


20 


APPENDIX  I 

CONTINUOUS  WAVE  LASER, 
CYLINDRICAL  BEAM  PROFILE 
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APPENDIX  II 


REPETITIVELY  PULSED  LASER, 
GAUSSIAN  BEAM  PROFILE 


PROGRAM  MPULSE ( INPUT * OUTPUT) 
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